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Abstract 

In this paper we use the method of matched asymptotic expan- 
sions in order to obtain a geometric motion as the singular limit of a 
nonlinear fourth order inhomogeneous equation. 
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1 Introduction 

Geometric evolution equations are often studied as the singular limit of 
reaction-diffusion equations. A typical example is the mean curvature flow 
equation 

V = H, (1) 

where H is the curvature of the front and V is the normal velocity to the 
front. It is obtained as the singular limit of the Allen-Cahn equation (Chen, 
1992) 

u t = Au ?ipu(u), (2) 

where e > is a small parameter and ip u is the derivative of the double-well 
potential ip(u) = ^(1 — u 2 ) 2 . If e —> 0, the solution u tends to the minimal 
values of ip, i. e. —1 and 1. Such a geometric evolution equation describes 
the motion of the boundary between regions where the limit equals —1 and 
1. In other words, the zero- level set of u evolves according to (pQ) when 
the curvature flow has a smooth solution. Mean curvature flow appears in 
the description of interface evolutions in a variety of physical models; this 
is related to the property that such a flow is the gradient flow of the area 
functional and therefore occurs naturally in problems where a surface energy 
is minimized. As a matter of fact, equation ([2]) is the gradient flow, suitably 



1 



rescaled, of the functional 



e IV7 |2 

2 |W| 



dx, 



where is a bounded open subset of R . In Loreti and March (2000), the 
authors are interested in the following functional 



^ |Vu| 2 + -ip(u) 



1 

dx H — 
£ Jn 



-eAu + -ip u (u) 



dx, (3) 



which is related to a conjecture that De Giorgi made about the approxima- 
tion, in the sense of T-convergence, of the Willmore functional, a functional 
depending on the curvature. They study such a functional in the field of 
image segmentation, that is, the recovery of smooth boundaries in a picture. 
The authors show that, in the case N = 3, the zero level set of the gradi- 
ent flow of Q, suitably rescaled, approximates a front rcfl that evolves 
according to the geometric motion 



V = H — 2A S H + AH ( K 



H- 



(4) 



where V is the normal velocity of the interface, A S H is the surface Laplacian, 
H and K are, respectively, the mean curvature and the Gaussian curvature 
of r at p(x), projection of x on T (i e. H = k± + hi and K = k\ ■ &2, 
being k\ and fc2 the principal curvatures of V). We observe that the term 
A S H is related to the Cahn-Hilliard equation (Cahn, Elliott and Novick- 
Cohen, 1996), since the zero level set of such an equation approximates, as 
e — >• and suitably rescaled, a front evolving with normal velocity propor- 
tional to the surface Laplacian of the mean curvature. 

Due to the increasing interest in reaction-diffusion equations with spatially 
inhomogeneous reaction terms deriving from Nakamura, Matano, Hilhorst 
and Schatzle (1999) and Lou (2003), this paper is concerned with the asymp- 
totic behaviour, as e — > 0, of the solution of the fourth order equation 

12 4 2 

ut = — w it b 2 (x)ip uu (u)w H — a(x)Va(x)Vw H — a 2 (x)Aw, (5) 



where 



w = b 2 l (x)ip u (u) — 2e 2 a(x)Va(x)Vu — e 2 a 2 {x)Au. 
It is the gradient flow of the following functional 



(6) 



n 



° 2 (X)! I v «| 2 + b 2 (x)-ip(u) 



dx 



(7) 



1 

+ - 

e Jn 



-ea (x)Au — 2ea{x)V a{x)V u + b (x)-ip u (u) 



dx, 
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where a{x) and b{x) are smooth and for some ao > 0, a{x) > ao and b(x) > 
Vi £ In the same line of Loreti and March (2000) we show, by using 
the formal method of matched asymptotic expansions, that the zero level 
set of a solution of ©, © approximates, as e — > and rescaling the time 
t — >■ et, a front evolving according to the geometric motion 



V = a 2 (x)H + 4Ha i (x) K 



+ 2a 3 (x)H 2 Va(x)Vd + 



a 5 (x) 
b(x) 



El 
4 

HA 



a(x) 



15 



Tr 2 a b (x)H 



V 



6(g) 
a{x) 



a 2 (x)A + 2a 2 (x)B + a(x)b(x)C I — ?T 

15 



1 



15 



a(x)b{x)D, (8) 



in the case iV = 3, where Vd and A,B,C,D are specified in sections O HI 
The paper is organized as follows. In section [21 we briefly describe the idea 
of the method of matched asymptotic expansions. Section [3] is devoted to 
the introduction of the problem and the instruments that we use to solve it. 
Finally, Section [J] contains the main calculations for the formal derivation 
of (ED. 



2 The method of matched asymptotic expansions 

In this section, we introduce the idea of the formal tool that we will adopt 
throughout the paper. Usually, the method of matched asymptotic expan- 
sions is used to solve singular perturbations problems. The general method 
with perturbation problems is to seek an expansion with respect to the 
asymptotic sequence {l,e,e 2 ,...},ase— > 

u e ~ Uq + eUi +e 2 U 2 + ..., e^O, 

for functions Uq, U\, U2, ■ ■ ■ which have to be determined. In a singular per- 
turbation problem, the regular methods produce an expansion that fails, at 
some point, to be valid over the whole domain. Then, the idea of the method 
is the following. The domain can be divided into two subdomains. On one 
subdomain, a solution is accurately approximated by an asymptotic series, 
obtained by treating the problem as a regular perturbation; we refer to it 
as the outer solution. On the other subdomain, the approximation cannot 
be accurate, since the perturbation terms in the problem are not negligible 
there. An approximation in the form of an asymptotic series is obtained 
there by handling this part of the domain as a separate perturbation prob- 
lem; we call it the inner solution. When the validity regions of the outer 
and inner expansions overlap, the outer and inner solutions are combined 
through a certain process (matching) in such a way that an approximate 
solution for the complete domain is found. 
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3 Formulation of the problem and preliminaries 



The gradient flow of the functional ([7]) is given by equations ©. We 
rescale the time t\ = et. The boundary conditions are obtained by taking 
the first variation of the functional ([7]) and provide that the first variation 
does not contain boundary terms. The problem to solve is then 
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-e 2 w - 2b 2 (x)ip uu (u)w + 4e 2 a(x)Va(x)Vtt> 
+2e 2 a 2 (x)Aw, in Ox(0,T) 

w = b 2 (x)Tp u (u) - 2e 2 a(x)Va(x)Vu 
-e 2 a 2 (x)Au, in Ox(0,T) 



(9) 



w 



n 



= 0, on d£l 

[ea 2 (x)Vu + §a 2 (x)Vu>] = 0, 



on 



on 



n being the outward unit normal to d£l. We denote the solutions of ([5]) 
and ([6]) with u and w respectively, keeping in mind their dependance on e. 
Moreover, we indicate the interface as 

r = u t > (r t x {t}) , 

where Tt = {x 6 R N : u(x,t) = 0}. We assume that the zero level set of a 
solution u of ([5]) divides the domain fi into two subdomains, Cl + and f2_. 
Following Cahn et al. (1996), Pego (1989) and Loreti and March (2000), we 
assume that u and w have the expansions 



and 



u(x, ti) ~ u° + eu 1 + e 2 u 2 + 



w(x, ti) ~ w° + ew 1 + e 2 ^ 2 + . . . 



(10) 



;n) 



away from the interface T. We refer to ()10p and (jlip as the outer expansions. 
The construction of an inner solution for x near T is more complicated: 
indeed, a correct scale must be decided in order to get inner expansions which 
match well the outer expansions. For this reason, in a small neighborhood 
of r, a stretched normal distance to the front is defined 



d(x,tx) 



(12) 



being d(x,ti) the signed distance from the point x £ f2 to T tl , such that 
d > for x S fi+ and d < for x € 0_. An ansazt is introduced as follows 
for the inner expansions 



u(x, h) = U(z, x, h) ~U° + eU 1 + e 2 U 2 + ... 



(13) 
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w(x, h) = W(z, x, ti) ~W° + eW l + e 2 W 2 + .... (14) 

As in Pego (1989) and Loreti and March (2000), it is required that the 
quantities depending on (z,x,ti) do not change when x varies normal to T, 
keeping z fixed. Then, we have the following setting 

V = d tl d{x,ti), 

being V the normal velocity of T in the t\ timescale, positive when the front 
moves toward Furthermore, indicating with Vc£(x,ii) the unit normal 
to r pointing toward 0+, we have 

Vm = -VdU z + V X U 

£ 

= V X U° + eV.f/ 1 + £ 2 V X U 2 + -Vd(U° + eU] + e 2 U 2 ) + ..., 

Au = A X U + -U z Ad+ \u zz 

= ^U Q zz + - £ {U° z Ad+Ul z )..., 

d tl u = d tl u + \vu z ; 

analogous equations are true for w and W. From Gilbarg and Trudinger 
(1977), in the case N = 3 

at/ , \ _ h(p(x)) k 2 {p(x)) 

{ ' l) l + kMx^dix,^) + i + k 2 ( P (x))d{x,t 1 y 

with p(x) projection of x on T. However, from ()12p and expanding on powers 
of e, it is easy to get 

Ad(x, t x )=H- £z(H 2 - 2K) + 0(e 2 ), (15) 

where H and K denote respectively the mean curvature and the gaussian 
curvature of the interface at p{x). From the definition of zero level set of T, 
the normalization conditions hold 

C/ fc (0,x,ti) = 0, A; = 0,1,.... (16) 

Finally, using the notation u±(x,ti) = lim r _ >0 ± u l (x + rV d(x , t\) , t\) , the 
following matching conditions are required: 

u\(xM)= hm U°(z,x,h), (17) 

z— >±oo 

lim (u± + zVdVu±)(x,h) = lim \J x (z, x, h), (18) 

z— >±oo z— s-±oo 

lim (u 2 ± + zVdVu\+\z 2 D 2 uV)(x,t 1 )= lim U 2 (z, x, tA (19) 

z— >±cxo 2 z— s-±oo 

where D denotes the directional derivative along Vd(x, t\) (see Loreti and 
March (2000) for further details). The expansions of w and W must satisfy 
analogous matching conditions. 
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4 Asymptotics 

In this section, we present a formal derivation of the interface equation (JSj) . 



4.1 The outer solution 

At the beginning, we start with the solution far from T. We notice that the 
minimizers of ([7D are expected to satisfy u(x) = ±1 almost everywhere; for 
this reason, we seek a solution u of ([5|) which is equal to —1 and 1 in 
and f2+ respectively, in the timescale t\ = et, as e approaching zero. The 
boundary conditions on <9f2 yield 

w * = 0, i = 0, 1,..., 
n ■ a 2 (x)Vw° = n ■ a 2 (x)Vw l = 0, 
n- (a 2 (x)Vu' ; + 2a 2 (x)Vw i+2 ) = 0, % = 0, 1, . . . . 



The substitution of the outer expansions (|10p and (jlip into ([9]) leads at zero 
order in e to the following problem 



(20) 



(21) 



f b 2 (x)i; uu (u )w = 

w° = 0, on dQ, 
n ■ a 2 (x)Vw° = 0, on <9S1 

A solution of problem (|20p . in accordance with our expectation, is 

^°(x,ti) = | +J' , u;°(x,t 1 ) = in 

Moreover, at first order in e, we have 

& 2 (X) [^uu{u°)w l +^uuu(u°)u 1 W°] = 

= 6 2 (x)VW n °) ul 
= 0, on dn 
n ■ o 2 (x)Vw 1 = 0, on <9f2 

whose solution is given by 

u 1 (x,t 1 ) = 0, w 1 (x,t 1 ) = 0. 

Going further in the analysis of second and third order, we have that u 2 and 
u 3 are null, hence the solution in the outer region is: 

u(x,h) = ±l + 0(e 4 ). 
4.2 The inner solution 

In the following, we deal with the solution near T. We substitute expansions 
(fT3j) and ([HI) into the equations of problem ([9]). 
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Zero order 



As already done in Cahn et al. (1996) and Loreti and March (2000), we 
assume W°(z,x,ti) = 0. At zero order in e the equation ([6]) becomes: 



a 2 {x)U° zz -b 2 {x)i(; u (U ) = 0. 



(22) 



The unique monotone increasing solution of (j22[) in accordance with the 
normalization condition (|16p and the matching condition (|17p . is 



U°(z, x,t\) = tanh 



6(x) 
a(x) 



Remark 4.1. We notice that the solution U° satisfies 



ir 



a(x) 



(23) 



(24) 



First order 

At first order in e, we have 

a 2 (x)W z 1 z - b 2 (x)^ uu (U Q )W l = 0, 

W 1 = ^(x^uuiU^U 1 - a 2 {x)U° z H - a 2 {x)U\ z . 
Let us define the operator L: 

LW 1 = a 2 (x)W z 1 z - b 2 {x)ij uu {U )W l . 

Clearly, from (|22|). we obtain 

Lf7 2 ° = 0, 

therefore we can assume 

~b(x) 



(25) 
(26) 



W 1 (z,x,t 1 ) = fj,(x,ti)U° 



a(x) 



with fi(x, t\) a bounded function to be determined. In this setting, equation 
(f26j) becomes 

LU 1 = - [a 2 {x)H + fi(x,h)] U°. (27) 
Multiplying by U z and integrating by parts, we obtain 



LU x U2dz 



LU^dz = 0. 



Thus 



LU L U^dz = -(a 2 (x)H(p(x),t 1 )) + n(xM)) / (U z ) 2 dz = 0. 
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Since J^^(U z ) 2 dz / 0, it follows [i(x,t\) = —a 2 (x)H(p(x),ti). Hence, 



W 1 (z,x,t l ) = -a 2 (x)H(p(x),t 1 )U^ 



6(a) 
a(x) 



From (J27D, we have LC/ 1 = 0, so that U^x,^) = c(x,ti)U° 
to (fl~6j) and to (fTBj) . the function c(x,ii) is null and therefore 

E^foz.ti) = 0. 



Kg) 

a(x) ' 



. Due 



Second order 

At second order in e, using the definition of the operator L, we deal with 
the equation 



LW 2 = -a 2 (x)W z i H -2a{x)Va{x)W z 1 Vd 
= -[a 2 (x)H + 2a(x)Va(x)Vd]W z 1 . 



(28) 



From (El), we have 



which gives 



b 2 {x)^ uu {U°)U 2 - 2a{x)Va(x)U° z Vd 
a 2 {x)A x U° + a 2 {x)zU° z {H 2 - 2K) - a 2 (x)U : 



2 

zzi 



LU 2 



-W 2 - 2a{x)Va{x)VdUl 



- a 2 (x)A x U° + a 2 (x)zU°(H 2 -2K). 
Calculations give the following expression for the laplacian of U° 



A x U l 



a(x) 



b{x) 
a(x) 



zU° + 



a 2 (x) 
¥{x) 



V 



Hx)_ 

a(x) 



Substituting it into (|29|) . we obtain 



-W 2 - 2a(x)Va(x)VdU° z 

a(x) 



+ a 2 {x)zU° z ^H 2 -2K 



b(x) 



A 



b(x) 
a(x) 



a (x) 
¥{x) 



b(x) 
a(x) 



z 2 U" zz . 



(29) 



(30) 



In order to obtain the function W 2 , we remind the matching condition (|19p 



lim W 2 (z,x,h) = 0. 

z— »±oo 
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The solution of (EEL is then 



W 2 = tf 2 (p(x),ti) \ a 2 {x)- + a(x,*i)J t/ 2 ° + a(x)H{p(x), h)V a(x)V dzU® , 



where a(x,ti) is a bounded function. Therefore 
a 2 (x)H 2 



2a 2 (x)K - a{x)HVa{x)Vd 



a 3 (x) 
b(x) 



A 



b(x) 
a(x) 



zU" 



[aH 2 + 2a{x)Va{x)Vd] £/° 



a 2 (x) 



a 2 (x) 
IPix) 



V 



6(g) 
a(x) 



z 2 C/° 



We remind that U 2 has to satisfy the matching condition (|19p 

lim C/ 2 (z,x,ti) = 0. 

z— >±oo 

With the same notation of Loreti and March (2000), if we set 

U 2 (z,x,t 1 ) = f(z,x,t 1 )U° z , 

we find that 



fz 
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a 2 {x){U®) 2 ' 



where g is such that 
a 2 (x)H 2 

9z - ' 



2a 2 (x)K - a(x)HVa(x)Vd 



a?(x) 
b(x) 



A 



b(x) 
a(x) 



+ 2 



[aH 2 + 2a(x)Va(x)Vd] (c/°) 2 
6(x) 



.a 3 (x) 



V 



a(x) 



z 2 U°(U° z ) 2 . 



Third order 

At third order in e, the first equation of problem Q becomes 



1, 



-VU" 



1 



(31) 



0\2 



W 1 + LW 3 - b 2 {x)i> uuu {U Q )U 2 W l 
a 2 (x)zWl{H 2 - 2K) + a 2 (x)W 2 H + a 2 (x)A x W 1 . 



(32) 



From now on, we omit to write the dependance on x of a and b. Here we 
need the expansion of the term A x W l . It may be obtained calculating 



AxW 1 = diviVW 1 ) = -div -V{abH)U° z 



+ 2zdiv 



a z V ( - ) HU°U° 
a 



9 



Using (|24p and collecting the different terms that we obtain from the calcu- 
lations, we can write 



^W 1 = AzU Q zz + BU% + Cz 2 (C/°) 2 + Dz 2 U°U° zz , 
where (from (Cahn et al., 1996), A X H = A S H) 

2 2 

A: = 2{Va) 2 H -2%r(Vb) 2 H -2^-VbVH + 2aVaVH 
b z b 

- ^—AbH + aAaH + 2^-VaSJbH - 2^-^-H, 
b b b 



2 

B : = -aHAa - 2aVaVH - 2^-VaVbH - %-HAb 

b b 

a 2 

- 2—VbVH - a 2 A s H, 
b 

(Vb) 2 b 
C := 2a^—--H - 4VaVbH - 2-(Va) 2 H, 
b a 

D := 2^(S/b) 2 H - 2VaVbH + 2-(Va) 2 H. 
b a 



Thus 



\VU° Z = ^a 2 HU° z + LW 3 + a 2 b 2 H^ uuu (XJ Q )U 2 XJ* 
+ a 4 HzU° z {H 2 - 2K) + a 2 W 2 H 

+ a 2 AzU zz + a 2 BU Q z +a 2 Cz 2 {U Q z ) 2 + a 2 Dz 2 U°U Q zz . (33) 



Taking into account the fact that ip uuu {u) = 12u, we multiply equation ([33 
by U z and integrate in z 

-I Z' + OO -I f + OO /•-* OC 

-V / (U°) 2 dz = -a 2 H / {U° z fdz+ / LW 3 U° z dz 

2 J -oo 2 J -CO J -oo 

/+oo 
U°U 2 (U°) 2 dz 
-oo 

/+oo 
zU^U^dz 
-oo 

r+oo r+oo 

+ a 2 H / W 2 U° y dz + a 2 A / zUWjz 



z~z — ' I ~" Z " zz" 

— oo J — oo 

/+oo r+oo 
(U° z ) 2 dz + a 2 C z 2 (U° z fdz 
-oo J — oo 

/+oo 
z 2 U°U^U^ z dz. 
-oo 
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Integration by parts yields 



We set 



therefore 



/_ 



-oo 

,3 TT 



LW 6 U"dz = 0. 



/+oo 
(U° z ) 2 dz, 
-oo 



i 1 r+oo 

±Vi 2 = -a 2 Hi 2 + 12a 2 b 2 H / U°U 2 (U°) 2 dz (34) 

/+oo 
zU^U^dz 
-oo 



-oo 

+ a 2 H I U° z W 2 dz 



oo 

/+oo 
zU° z U° zz dz + a 2 5i 2 
-oo 
-oo 

+ a 2 C / z 2 «7°) 3 cfe 



/+oo 
z 2 (U^f 
-oo 

/+oo 
z 2 C/°?7 2 ?7 2 2 ^. (35) 
-oo 



We calculate explicitly the integrals in the above expression. From integra- 
tion by parts, it is easy to obtain 



+°° ?'n 

zU° z U° zz dz = 

Z ZZ O 

oo z 



and consequently 



/+oo -I /-+00 -I 

zU°(U° z ) 2 dz = --j / zU° z U° zz dz = - A 



a 

—i2- 
46 
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Afterwards, due to (JMD and ([31 

r+oo 



U°U 2 (U° z ) 2 dz 



+ 



_la 
~2b 
1 1 
6 ab 
1 1 
6 ab 
1 1 
6 aft 
1 1 
6 ab 



— oo 
+00 



U° z J{U° z fdz 
gUUz 



oo 

— oo 
2 tlt2 



2a 2 iT - aHVaVd 



-A 



zU Q {U° z fdz 



1 1 2 a 3 
6 a& 6 



/+oo 
U°(U°) 2 dz 
-oo 

Hoo 

z 2 ([/°) 2 ([/°) 2 ^. 



V 



The integral J_ °° U°(U z ) 2 dz is null since the integrand is an odd function 
with respect to z, while 



/ 



+oo 2 3 

z\U^{U«Tdz=*- a - 
oo 45 o 



Thus 

P+OO 



U°U 2 (U° z ) 2 dz 



J_ ^2 

lo^vr 2 
3 b 3 45 



a 2 if 2 



2a 2 iT - aHVaVd - ^- A ( - 



1 2 



Furthermore 

P + OO 



= / + °° (a 2 + §VaVd) + zU° zz ) 



-oo 
r+oo 



f/2 



+ 



aH 2 U° z U° zz dz 



-oo 
r2 



+ ^ (a 2 + fvaVd) / 



°) 2 



oo 
+oo 



U z U° zz dz 



a 2 H 2 aHX7aX7d\ 



+ 
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Finally 



— oo 



+00 



and 



—00 



z 2 U°U° z U° zz dz 



45 



The substitution of the above integrals into ([8]) gives the following interface 
equation 



if N = 3. 

Conclusions 

In this paper, we have studied a geometric law considering the gradient flow 
of the functional ([7]) , built following the variational motivation used in Loreti 
and March (2000). We have proved, by means of formal asymptotics, that 
such a motion may be approximated in a suitable sense by an inhomoge- 
neous fourth order parabolic equation. It involves the mean and Gaussian 
curvatures and the surface Laplacian of the mean curvature of the evolving 
interface. As in Nakamura et al. (1999), the motion equation arising from 
([5]) and involves drift terms, despite the absence of drifts in the original 
equation. 
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